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We consider a simple linear reversible isomerization reaction A ^ B under subdiffusion described 
by continuous time random walks (CTRW). The reactants' transformations take place independently 
on the motion and are described by constant rates. We show that the form of the ensuing system 
of mesoscopic reaction-subdiffusion is somewhat unusual: the equation giving the time derivative of 
one reactant concentration, say A(x,t), contains the terms depending not only on AA, but also on 
A_B, i.e. depends also on the transport operator of another reactant. Physically this is due to the 
fact that several transitions from A to B and back may take place at one site before the particle 
jumps. 

PACS numbers: 05.40.Fb, 82.40.-g 
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There are several reasons to discuss in detail the struc- 
ture of mesoscopic kinetic equations describing the be- 
havior of a simple reversible isomerization reaction A ^ 
B under subdiffusion. 

Many phenomena in systems out of equilibrium can 
be described within a framework of reaction-diffusion 
equations. Examples can be found in various disciplines 
ranging from chemistry and physics to biology. Both 
reaction-diffusion systems with normal and anomalous 
diffusion have been extensively studied over the past 
decades. However, for the latter, a general theoretical 
framework which would hold for all kinds of reactions 
is still absent. The reasons for subdiffusion and there- 
fore its properties can be different in systems of different 
kind; we concentrate here on the situations when such 
subdiffusion can be adequately described by continuous- 
time random walks (CTRW). In CTRW the overall par- 
ticle's motion can be considered as a sequence of jumps 
interrupted by waiting times, the case pertinent to many 
systems where the transport is slowed down by obsta- 
cles and binding sites. In the case of anomalous diffu- 
sion these times are distributed according to a power law 
lacking the mean. The case of exponential distribution, 
on the other hand, corresponds to a normal diffusion. 
On the microscopic level of particles' encounter the con- 
sideration of subdiffusion does not seem to be problem- 
atic, although it has posed several interesting questions 
[E S IS BL 01 • However, these microscopic approaches 
cannot be immediately adopted for description of spa- 
tially inhomogeneous systems, which, in the case of nor- 
mal diffusion, are successfully described within the frame- 
work of reaction-diffusion equations. To discuss such be- 
havior under subdiffusion many authors used the kind 
of description where the customary reaction term was 
added to a subdiffusion equation for concentrations to 
describe such phenomena as a reaction front propagation 



or Turing instability @, 0, & i, E3, EH • 

The results of these works were jeopardized after it was 
shown in Ref. [12| that these procedure does not lead to a 
correct description even of a simple irreversible isomeriza- 
tion reaction A — * B. The transport operator describing 
the subdiffusion is explicitly dependent on the properties 
of reaction, which stems from an essentially very simple 
observation that only those particles jump (as A) which 
survive (as A). 

The properties of the reaction depend strongly on 
whether the reaction takes place only with the step of 
the particle, or independently on the particles' steps, 
and moreover, whether the newborn particle retains the 
rest of it previous waiting time or is assigned a new one 
[Hil H3 |. Here we consider in detail the following sit- 
uation: The A ^ B transformations take place inde- 
pendently on the particles' jumps; the waiting time of 
a particle on a site is not changed by the reaction, both 
for the the forward and for the backward transformation. 
As a motivation for such a scheme we can consider the 
reaction as taking place in an aqueous solution which 
soaks a porous medium (say a sponge or some geophys- 
ical formation). If sojourn times in each pore are dis- 
tributed according to the power law, the diffusion on the 
larger scales is anomalous; on the other hand, the reac- 
tion within each pore follows usual kinetics. We start 
by putting a droplet containing, say, only A particles 
somewhere within the system and follow the spread and 
reaction by measuring the local A and B concentrations. 

Stoichiometry of the chemical reaction implies the ex- 
istence of a conservation law. In the case of the A ^ B 
it is evident that the overall number of particles is con- 
served. If the isomerization takes place independently on 
the particle's motion, then the evolution of the overall 
concentration C(x,t) — A(x,t) + B(x,t), where A(x,t) 
and B(x,t) are the local concentrations of A and B par- 
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tides respectively, is not influenced by the reaction, and 
has to follow the simple subdiffusion equation 

C(x,t)=K a0 D 1 t - a AC 

as it should be. On the other hand, neither the result 
of the treatment in Ref.[l5[ nor the result of Ref.pjJ 
reproduce this behavior which is a consequence of the 
fundamental stoichiometry. In the work [151 ] (where two 
of the authors of the present report were involved) it 
was implicitly assumed that the back reaction can only 
take place on a step of a particle, without discussing this 
assumption. The corresponding description lead to the 
expressions which could not be cast in a form resembling 
the reaction-diffusion equations at all. The more general 
approach of Ref . (l6l | , definitely correct for irreversible re- 
actions, also fails to reproduce this local conservation law 
and thus is inappropriate for the description of reversible 
reactions under the conditions discussed. According to 
Ref. [l4| the approach of Ref. 16] implies that the waiting 
time after each reaction is assigned anew, which makes a 
large difference in the reversible case. 

As a step on the way to understanding the possible 
form of the reaction-subdiffusion equations we consider in 
what follows the simplest linear reversible scheme where 
each step can be explicitly checked. We show that the 
form of the corresponding equations is somewhat un- 
usual, which emphasizes the role of coupling between the 
reaction and transport in reaction-subdiffusion kinetics. 
Actually, the equation giving the time derivative of one 



reactant concentration, say A(x,t), contains the terms 
depending not only on AA, but also on AB, i.e. depends 
also on the transport operator of another reactant. Phys- 
ically this is due to the fact that several transitions from 
A to B and back may take place at one site before the 
particle jumps. This dependence disappears only in the 
Markovian case due to vanishing of the corresponding 
prefactor. 

Following the approach of Ref. [H, [H[ we describe the 
behavior of concentrations in the discrete scheme by the 
following equations: 

Ai{t) = -h(t) + i/<_i(t) + \h+i{t) - k x A{t) + k 2 B{t) 
Bi{t) = -Ji{t) + i J<_i(t) + \j l+ i{t) + k x A{t) - k 2 B{t) 

where I{t) is the loss flux of A-particles on site i and J(t) 
is the corresponding loss flux for B-particles at site i. In 
the continuous limit the equations read as 

a 2 

A(x,t) = —AI{x,t)-kxA{x,t) + k 2 B{x,t) (1) 
a 2 

B(x,t) = —AJ(x,t) + k x A{x,t)-k 2 B{x,t). (2) 

We now use the conservation laws for A and B particles 
to obtain the equations for the corresponding fluxes. The 
equations for the particles' fluxes on a given site in time 
domain (the index i or the coordinate x is omitted) are: 



/(*) - i>(t)P AA {t)A(0) + / V(i - t')P AA {t - t') I{t')+k 1 A{t')-k 2 B{t')+A[t') 



iP(t)P BA (t)B{0) + i>{t - t')P BA {t - t') J(t') - ktAit') + k 2 B(t') + B(t') 
Jo 



for A-particles, and 



(3) 



J(t) = ip(t)P BB (t)B(0) + / i>(t-t')P BB {t-t') J(t')-kiA(t') + k 2 B(t')+B{t') 



i/>(t)P AB (t)A(0) + / tl>(t-t')P AB (t-t') J(t') + fciA(t') - k 2 B(t') + A{t') 
Jo 



for B-particles. 



The explanation of the form of e.g. Eq.([3]) is as fol- 
lows: An A-particle which jumps from a give site at time 
t either was there as A from the very beginning, and 
jumps as A probably having changed its nature several 
time in between, or came later as A and jumps as A, or 
was there from the very beginning as B and leave the 
site as A, etc. Here P AA , P AB , P BA and P BB are the 
survival/transformation probabilities, i.e. the probabil- 
ity that a particle coming to a site as A at t = leaves 



it at time i as A (probably having changed its nature 
from A to B and back in between), the probability that 
a particle coming to a site as A at t = leaves it at time 
t as B, the probability that a particle coming to a site as 
B at t — leaves it at time t as A, and the probability 
that a particle coming to a site as B at t = leaves it at 
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time t as B: 

P AA {t) 
P BA {t) 
P BB {t) 
Pab it) 



I + fc l . e -(kx+k 2 )t 



ki + k 2 ki + k 2 

k 2 k 2 
ki + k 2 ki + k 2 
k 



-(k 1 +k 2 )t 



fci + k 2 ki + k 2 

k\ k\ 
k\ + k 2 ki + k 2 



-(k 1 +k 2 )t 



(4) 



These are given by the solution of the classical reaction 
kinetic equations 



A(t) = -kiA(t) + k 2 B(t) 
B(t) = kiA(t) - k 2 B(t). 



(•5) 



The values of P AAl Pab are given by the solutions 
P AA (t) = A(t) and P A s{t) = B(t) under initial condi- 
tions A(0) = 1, B(0) = 0, and the values of Pba and 
Pbb are given by P B A{t) = A(t) and P B B{t) = B(t) 
under initial conditions A(0) = 0, B(0) = 1. 
In the Laplace domain we get: 

I(u) = %l) X (u)[I{u) + kiA(u)-k 2 B{u)+uA(u)] 

+ip 2 {u) [J(u) - kiA(u) + k 2 B(u) + uB(u)} 

J[u) = i/) 3 (u) [J(u) - k x A(u) + k 2 B(u) + uB(u)] (6) 
+ip4(u) [I(u) + kiA(u) - k 2 B(u) + uA{u)] 

where ipi(u),ip 2 (u),ip3(u) and ip4,(u) are the Laplace 
transforms of rp^t) = ip(t)P AA (t), ip 2 (t) = ip(t)P BA (t), 
fait) = ip(t)P B B(t) and ^ 4 (t) = ip{t)P AB {t), respec- 
tively. 

Using shift theorem we can get the representations of 
ipi in the Laplace domain. They read: 



1pi{u) 



fci + fc 2 
fc 2 

fci + k 2 

fci 
fci + k 2 

fci 
fci + k 2 



ip{u) 
tp{u) 
ip(u) 
ip(u) 



fcl + fc - 2 

k 2 

fcl + fc - 2 

fc 2 
fcl + fc - 2 

fcl 
fcl + fc 2 



1p(u + kl +fc 2 ) 
ip(u + fcl + fc 2 ) 
tp(u + fcl + fc 2 ) 
ipiu + h +k 2 ) (7) 



The system of linear equations for the currents, Eqs.©, 
then has the solution 

I(u) = an(u)A(u) + ai 2 {u)B{u) 
J(u) = a 2 i(u)A(u) + a 22 {u)B(u) 

with the following values for the coefficients: 
1 1 



an 



fci + fc 2 1 + <jnp — ip — cf) 
x [-<ptp (fcifc 2 + u(ki + fc 2 ) + fci) 
+ (j>ki(u + fci + fc 2 ) + ipk 2 u] , 



0.21 



fcl 



fci + fc 2 1 + (pip — ip - 4> 

X [^(fci + fc 2 ) + (ip- 4>)u - (fcl + fc 2 ) 



and with the two other coefficients, ai 2 and a 22 differing 
from a 2 i and an by interchanging fci and fc 2 . Here ip = 
ip(u) and <p = ip(u + fci + fc 2 ). 

For the exponential waiting time density ip(t) = 
t~ 1 exp(— tjr) the corresponding values get 



an = a 22 = 1/r 
ai2 = a 21 = 0, 



and the system of equations for the concentrations in 
the continuous limit, Eqs.Q, reduces to the customary 
system of reaction-diffusion equations. For the case of 
the power-law distributions ip(t) — t~ l ~ a the Laplace 
transform of the waiting time PDF is ip(u) ~ 1 — cu a for 
small u, with c = t°T(1 — a), so that 



an 



a 22 



[fcaw 1 a + k 1 {u + k 1 +k 2 ) 



l-a 



a 2 i = 



fci + k 2 
cfci(fci + k 2 ) - cu{k\ + fc 2 )] 
r- 1 

[km 1 -* + k 2 (u + k 1 + k 2 ) 1 - a 
cu(ki + fc 2 )] 

[fciu 1 -" -ki(u+h + k 2 y- a 



fci + fc 2 

cfc 2 (fci + fc 2 



fcl + fc 2 

Cfci(fci +fc 2 )] 

r .-l 



a 12 



fci + fc 2 " 
+ cfc 2 (fci+fc 2 )] 



hu 1 -* - k 2 ( u + k! + k 2 y- a 



Now we turn to the case of long times and relatively slow 
reactions, so that all parameters, u, fci and fc 2 can be 
considered as small. In this case, for a < 1, the leading 
terms in all these parameters are the first two terms in 
each of the four equations, and the other terms can be 
neglected. In the time domain the operator correspond- 
ing to u 1_Q is one of the fractional derivative oD\~ a , 
and the operator corresponding to (u + fci + fc 2 ) 1_Q is 
the transport operator of Ref . [la] . "(fci + fc 2 ) with 
oT t 1 ~ a (k) = e~ ht oDl~ a e kt , Introducing the correspond- 
ing equations for the currents into the balance equations 
for the particle concentrations we get: 
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A(x,t) = K a 
+ K a 

B(x,t) = K a 
+ K a 



ki + k 2 

k 2 
h + k 2 

ki 
ki + k 2 

fci 
k\ + k 2 



2 7-\l — CK 

U t 



r>i-a 



[ ) 



ki + k 2 

k 2 
ki + k 2 

fcl 

ki + k 2 

fc 2 
k\ + k 2 



T*- a (kx + k 2 ) 

r t 1 - Q (fc 1 + fc 2 ) 



AA(x,t) 

AB(x, t) - kiA(x, t) + k 2 B(x, t) 
AA(x,t) 

AB(x, t) + k\A(x, t) - k 2 B(x, t). 



(8) 



(9) 



Note also that the equation for C(x, t) = A(x, t) + B(x,t) 
following from summing up the Eqs.([8]) and @ is a simple 
subdiffusion equation 

C(x,t) = K a D\- a AC 

as it should be. On the other hand, neither the result of 
the treatment in Ref . [l5[ nor the result of Ref . [l6[ repro- 
duce this behavior which is a consequence of the funda- 
mental conservation law prescribed by the stoichiometry 
of reaction. 

Note that this system still holds for a = 1 when both 
the fractional derivative QD t ~ a and the transport oper- 
ator T t 1 ~ a (ki + k 2 ) are unit operators. In this case the 
usual system of reaction-diffusion equations is restored: 

A(x,t) = KAA(x,t) - kxA{x,t) + k 2 B(x,t) 
B(x,t) = KAB(x, t) + kiA(x, t) - k 2 B(x, t). 

Let us summarize our findings. We considered the sys- 
tem of mesoscopic (reaction-subdiffusion) equations de- 
scribing the kinetics of a reversible isomerization A ^ B 
taking place in a subdiffusive medium. When the wait- 
ing times of the particles are not assigned anew after 
their transformations (i.e. when the overall concentra- 
tion of reactants is governed by the simple subdiffusion 



r 



equation), this reaction is described by a rather unusual 
system of reaction-subdiffusion equations having a form 
which was up to our knowledge not discussed before: 
Each of the equations, giving the temporal changes of 
the corresponding concentrations, depends on the Lapla- 
cians of both concentrations, A and B (not only on the 
same one, as in the case of normal diffusion). This is 
a rather unexpected situation especially taking into ac- 
count the fact that our reaction is practically decoupled 
from the transport of particles. The form reduces to a 
usual reaction-diffusion form for normal diffusion (due to 
cancellations). It is important to note that the physical 
reason of the appearance of such a form is the possibility 
of several transformations A — > B — > A — > B... during 
one waiting period, and that such possibilities have to 
be taken into account also for more complex reactions 
including reversible stages. 
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